Multifractality of eigenfunctions in spin chains 
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We investigate different one-dimensional quantum spin-i chain models and by combining ana- 
lytical and numerical calculations prove that their ground state wave functions in the natural spin 
basis are multifractals with, in general, non-trivial fractal dimensions. 
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- One- dimensional quantum spin chains are among the 
oldest and the most investigated fundamental models in 
physics. Introduced as toy- models of magnetism they 
quickly became a paradigm of quantum integrable models 
(see e.g. A prototypical example is the XYZ 

Heisenberg model [l[ for N spins 
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and its various specifications for different values of pa- 
rameters. <Jn ,v,z are the Pauli matrices at site n. 

In the natural basis of z-components of each spin, \a) = 
|<ti , . . . , <7/v) where Oj = ±1, any Hamiltonian of N spins- 
\ is represented byaMxM matrix with dimension M = 
2 N . Many different methods were developed to determine 
exact spectra of such matrices The calculation of 

eigenfunctions is more involved. A wave function of N 
spins- \ in the spin-z basis can be written as 



(2) 



where the summation is taken over all M — 2 config- 
urations with <jj = ±1. In general, coefficients "J g can 
be found only after the matrix diagonalization which for 
large N is a hard numerical problem. Even in integrable 
cases eigenfunctions of spin chains look erratic (cf. fig- 
ures below) and their structure is not well understood. 

The purpose of this letter is to prove that ground state 
(GS) wave functions for different one-dimensional spin 
chain models are multifractals in the spin-z basis. 

- Multifractality is a general notion introduced to char- 
acterize strong and irregular fluctuations of various quan- 
tities For eigenfunctions like in (0) one uses the 
following definition (see e.g. Q and references therein). 
Let Sn(q,M) be the Renyi entropy for an eigenfunction 
([2]) of a matrix of finite size M 
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Fractal dimensions, D q , are defined from the behaviour 
of the Renyi entropy ([3]) in the limit M — > oo [9| 



lim 
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The case when D q is a non-linear function of q corre- 
sponds to a multifractal irregular behaviour. 

Fractal dimensions give a concise description of wave 
function moments and are important characteristics of 
eigenfunctions but it seems that they were overlooked in 
previous studies. For example, the Renyi and Shannon 
entropies are calculated for GS wave functions of certain 
spin chains in [I3|-[l3| but terms linear in In M which de- 
termine D q in ^ were regularly ignored and only next- 
to-the leading terms have been investigated as it is usual 
in conformal field theories. To the best of authors' knowl- 
edge only a recent paper (l^ | briefly mentioned the exis- 
tence of multifractality in a spin chain model. 

The simplest method to find fractal dimensions is the 
direct numerical calculation of the GS wave function for 
different number of spins and a subsequent extrapolation 
of the Renyi entropy for large M. Definition (j4|) is well 
suited for positive q. For many problems (but not for all) 
fractal dimensions can be calculated also for negative q 
[IH; 16 1 . Of course, if certain coefficients in ^) are zero 



with normalized coefficients ^g, J2s l^<?| 2 = T 



due to an exact symmetry, they are not included in the 
calculation of the Renyi entropy for q < 0. 

In what follows we investigate various specifications of 
the Heisenberg model ([1} and by combining numerical 
and analytical methods demonstrate that the multifrac- 
tality of GS is a generic property of all of them. We 
choose 7 > and a > to ensure off-diagonal terms of 
Hamiltonian matrices ([T]) to be non-positive which, by 
the Perron-Frobenius theorem, implies that coefficients 
^3 in ((2} for the GS wave function are non-negative. 
Other parameters are such that GS wave functions of 
the XY models (with A = 0) are ferromagnetic and for 
the XYZ models they are anti-ferromagnetic. 

- The quantum Ising model in transverse field 17| is 
a standard model of quantum phase transitions [18(. It 
corresponds to the Hamiltonian (JXJ) with A = a = and 
7 = 1. Its spectrum can be found analytically by the 
Jordan- Wigner transformation Q and coefficients \^g\ 2 
are given by the determinant of N x N matrices Q , [ll| . 
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FIG. 1. (color online). Fractal dimensions of GS for the 
quantum Ising model. Red line: A = 1.6, black line: A = 1, 
blue line: A = 0.4. Inset: GS coefficients for A = 1 and 
N = 11. The abscissa axis here and in other figures is integer 
binary code for 8, x = ^ =i 2 n_2 (l + o n ). 



Fractal dimensions of GS wave function for this model 
computed numerically from linear extrapolation of the 
Renyi entropy with N = 3 — 11 are presented in Fig. [T] 
for a few values of transverse field A. The curves of D q as 
a function of q for all models have the same characteristic 
form as for other fractal measures [7|. In particular, when 
q — > ±oo they tend to well defined limits D± OQ . 

Generalizing results of ll|, one gets the exact expres- 
sions for limiting values D± OQ and for -D1/2 



D 1/2 (A) = !-£>«, (i 



A — cos u 



In 1 ± , du, 5) 

L v^R(A^)J 

R(X,u) = l-2Acosu + A 2 . 



These formulas prove that fractal dimensions of quan- 
tum Ising model are non-trivial. In Fig. [2] these exact 
expressions are plotted together with numerically calcu- 
lated points for different A. D± oc arc obtained by a ht 
D±oo + a/q + b/q 2 fo large q parts of curves similar to 
Fig.[TJ The good agreement between Eqs. and numer- 
ics shows that though we obtain fractal dimensions from 
relatively small number of spins our results are reliable. 

The above curves are qualitatively the same for non- 
critical and critical (that is A = 1) quantum Ising model. 
Nevertheless, as illustrated in the inset of Fig. [5] the sum 
of and -D_oo has clear singularity in the critical point 
in accordance with the relation 



I?_ 00 (A)+^ OC (A) 



2, |A|<1 
|A|>1 
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which follows from ([5]). It means that criticality can be 
observed in fractal dimensions of GS. 

- The XY model is a specification of (JXJ) with A = 
a = and 7 7^ 1. Similar to the quantum Ising model 
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FIG. 2. (color online). Exact fractal dimensions ([S]) for the 
quantum Ising model: (blue line), D1/2 (black line), and 
D-oo (red line). Results of numerical calculations are indi- 
cated by symbols of the same color. Inset: Doo + calcu- 
lated numerically (black triangles) in comparison with Eq. (J6]) 
(red solid line). 



this model is also integrable by the Jordan- Wigner trans- 
formation 



ll| but the structure of its GS is more 
complicated. An interesting special case is A — A f where 
A/ = \Jl — 7 2 . It is known [19| that at that field the XY 
model has two exact factorizcd GS wave functions 



TT(cos0|l> n ±sin0|-l> n ), cos 2 20=^— J -. (7) 
, J- + 7 



For states with definite parity and A = A/ fractal dimen- 
sions are described by a formula 



D q = - 



ln(cos 2 9 9 + sm 2q 6) 
(q- l)ln2 



(8) 




FIG. 3. (color online). Fractal dimensions of GS for the XY 
model with anisotropy 7 = 1.4. Red line: A = 1.6, black 
line: A = 1, blue line: A = 0.4. Dashed black line shows 
for comparison the exact fractal dimensions (JH]) for A = 0.8 
and 7 = 0.6. Inset: GS coefficients for A = 0.4, 7 = 1.4 and 
N = 11. 
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FIG. 4. (color online). Asymptotic fractal dimensions for the 
XY model with anisotropy 7 = 1.4 versus transverse field. 
Red and blue lines indicate D± in @. Black line is the con- 
tribution of the Neel configuration (|10[) . When contributions 
become sub-dominant they are indicated by dashed lines of 
the same color. Blue squares and red triangles are respectively 
D-oo and D x calculated as for the Ising model. 
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FIG. 5. (color online). Fractal dimensions for the quantum 
Ising model in transverse (A) and longitudinal (a) fields. Red 
line: a = 2, A = 2, black line: a = 1.4, A = 2, blue line: 
a = 0.6, A = 0.4. Inset: GS coefficients for a = 1.4, A = 2, 
and N = 11. 



indicated for 7 = 0.6 and A = 0.8 by dashed black line 
in Fig. [3] This example proves that at least at the fac- 
torizing field fractal dimensions of the GS wave function 
do exist and correspond to the well-investigated case of 
binomial measures [6J. 

As for A 2 + 7 2 < 1 there exist many crossings of low- 
est states with different parity, for numerical calculations 
(performed as in the Ising model) we choose A and 7 out- 
side the unit circle, A 2 +7 2 > 1. The results are presented 
in Fig. [3] and are qualitatively similar to the Ising model. 

One may argue that the limiting values, and i3_oo 
as in the quantum Ising model should correspond to con- 
figurations with, respectively, all spins up and all spins 
down, and, consequently, are expressed similar to (J5J) as 
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cos it) 2 + 7 2 sin 2 u. But for small 



with i?_(A,7,w) = (A 
A the minimal contribution is instead given by the anti- 
ferromagnetic Neel configuration with alternating spins, 
a n = (—IT™. Using the asymptotics of the block Toeplitz 
matrices [20| we get 

-D N eel(A, 7) = 



(10) 
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where R+(X, 7, u) = (A + cosu) 2 + 7 2 sin 2 u. This result 
for 7 = 1.4 is presented in Fig.|4]by the black line. When 
L>Ncci > £>_, .D-oo = L>Ncoi, otherwise £>_oo = D— For 
7 = 1.4 these curves intersect at A ~ 0.4982 and -D-oc 
has the form indicated in Fig. 2] by solid blue and black 
lines. Numerical results agree well with this prediction. 

- The Ising model in transverse and longitudinal fields 
is obtained by adding to the quantum Ising model a longi- 
tudinal field a. For non-zero a the only known integrable 



case corresponds to A = 1 [21|. This model attracts re- 



cently wide attention as certain consequences of its inte- 
grability have been checked experimentally in the cobalt 
niobate ferromagnet (22[. In Fig.[S]the fractal dimensions 
for a few values of both fields are presented. 

- The XXZ model in zero fields is a particular case 
of the Heisenberg model Q with 7 = A = a = and 
A 7^ 0. Due to the conservation of the z component 
of the total spin, S z — J2n a ni its Hamiltonian can be 
diagonalized in subspace with fixed S z . The model is 
soluble by the coordinate Bethe anzatz 0, [13, and 
has a rich phase diagram (see e.g. Q). 

As a reference we use A = — | called combinato- 
rial point. From the Razumov-Stroganov conjecture 
(25| proved in [26] it follows that at such A and odd 
N = 2R + 1 the following statements are valid: (i) the 
GS energy is —3^/4, (ii) the largest coefficient in the 
expansion (the one for the Neel configuration) equals 



3 fl / 2 2 ■ 5 . . . (3g - 1) 
2 R 1 -3...(2R- 1) 



(11) 



(iii) the smallest coefficient corresponding to a half con- 
secutive spins up and other spins down is ^^ in = 
^^.Ar, and (iv) J2s ^ = 3 R/2 - Here A R is the num- 
ber of alternating sign matrices [27j . 

These formulas prove that for A = — h fractal dimen- 
sions Doo and -D1/2 are explicitly known 
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When R -> 00, hiA R = ij 2 ln(3\/3/4) + O(R). Such 
quadratic in N behaviour is a particular case of the 
emptiness formation probability of a string of n aligned 
spins with n ~ N [28]. This asymptotics means that 
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FIG. 6. (color online). Fractal dimensions for the XXZ model 
in zero fields. Red line: A = —0.25, black line: combinatorial 
point A = —0.5, blue line: A = —1.6. Inset: GS coefficients 
for A = —0.5 and N = 13. Dashed line indicates theoretical 
prediction for Doo (|12|l at A = —0.5. Black circle is the value 
of -D1/2 l|12p at this point. 
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FIG. 7. (color online). Fractal dimensions for the XYZ model 
with anisotropy 7 = 0.6. Red line: combinatorial point A = 
-0.32, black line: A = -0.5, blue line: A = -1. Inset: GS 
coefficients for A = —0.5, 7 = 0.6, and N = 12. 



negative moments of the GS wave function in anti- 
ferromagnetic case require a scaling different from 
and hence will not be considered here. The fractal di- 
mensions for a few values of parameter A are presented 
in Fig. [6l The numerical calculations were performed by 
extrapolation of the Renyi entropy separatively for odd 
and even N = 3 — 13. With available precision, fractal 
dimensions for odd and even N are the same but sub- 
leading terms in the Renyi entropy ([3]) are different. 

- The XYZ model differs from the XXZ model by 
anisotropy 7 7^ in the (x, y) plane. In zero fields its 
GS wave function has been found in [H[ . A soluble exam- 
ple with factorized GS ([7])isa = 0,A = A/ where 



X f = v/(l-A) 2 -7 2 , cos 2 2g= 7 (13) 

1 + 7 — A 

At this field GS wave function corresponds to the bino- 
mial measure and fractal dimensions are given by ([5)1. 

The combinatorial point for the XYZ model in zero 
fields where an additional information about GS is known 
(or conjectured) is A = (j 2 - l)/2 29]. In Fig. 0we 
present fractal dimension in the zero fields XYZ model 
with 7 = 0.6 for a few values of A including the com- 
binatorial point. Qualitatively the curves are similar to 
the XXZ model and to D q with positive q for the XY 
models. 

- Summary. Wave functions are the most fundamen- 
tal objects of any quantum- mechanical model. For many- 
body problems their structure is complicated and numer- 
ous questions remain open. We consider here practically 
all standard one-dimensional spin-^ models and demon- 
strate that their GS wave functions in the natural spin-z 
basis are multifractals with, in general, non-trivial fractal 
dimensions. For special values of parameters and/or cer- 



tain dimensions we get exact analytical formulas which 
prove rigorously the existence of fractal dimensions. In 
other cases we rely on numerical calculations. The multi- 
fractality in spin chains is a very robust phenomenon. It 
exists for integrable and non-integrable models, for ferro 
and anti-ferro magnetic states, as well as for critical and 
non-critical systems. 

A common point of all these models is that their MxM 
Hamiltonian matrix in spin-z basis is such that in each 
row and column there exist only K ~ In M non-zero ma- 
trix elements of the same order. As M 3> K this specific 
form resembles a tree structure with branching number 
K. As K — > 00 when M — > 00 and it is known [3(j, [3l| 
that on a tree with large branching number (with other 
parameters fixed) the Anderson localization is unlikely, 
states are delocalized. On the other hand, the full ergod- 
icity on a tree is, in general, also improbable [3X1 ] . The 
remaining possibility corresponds to delocalized but not 
ergodic states, i.e. to multifractality, which may explain 
its ubiquity. The multifractality is related finally not 
with explicit randomness but with internal complexity of 
models considered. Such arguments are not restricted to 
one dimension and/or spin chains but can be applied to 
various many-body problems with local interactions and 
we conjecture that multifractality (i.e. a non-trivial scal- 
ing of different quantities with the number of particles) is 
a generic property of a large class of many-body models. 
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